Symmetry Violation in a Superconducting Film 
with a Square Array of Ferromagnetic Dots. 
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We study the equilibrium state of a superconducting film covered with a regular square array 
of perpendicularly magnetized magnetic dots. The dots induce vortices in the film directly under 
them and antivortices between the dots. We show that the symmetry of the dot array is sponta- 
neously violated by the vortices. The positions of the vortices and the antivortices depend on the 
magnetization and the size of the dots. 



Heterogeneous structures made of ferromagnetic (FM) 
and superconducting (SC) pieces have attracted a great 
deal of interest both experimentally [1-7] and theoreti- 
cally in recent years [8-15]. In contrast to conventional 
systems, these structures allow the coexistence of ferro- 
magnetism and superconductivity, if FM and SC order 
parameters are separated in space by a thin insulating 
layer [8]. Strong interaction of the FM and SC systems 
not only gives rise to a new class of novel phenomena and 
physical effects, but also shows the important technolog- 
ical promise of devices whose transport properties can be 
easily tuned by comparatively weak magnetic fields. 

Several realizations of such systems were proposed. 
However, so far only sub-micron magnetic dots covered 
by thin SC films have been prepared and studied [1—7]. 
The experimental samples of FM-SC hybrid systems were 
prepared by means of electron beam lithography and lift- 
off techniques [17]. Both in-plane and out-of-plane mag- 
netization was experimentally realized. The dots with 
magnetization parallel to the plane were fabricated from 
Co, Ni, Fe, Gd-Co and Sm-Co alloys. For the dots with 
magnetization perpendicular to the plane, experimenters 
used Co/Pt multilayers. The FM dots were deposited on 
thin SC films made of either Nb or Pb, whose transition 
temperatures are around 7— 10 K. In these experiments, 
the effects of commensurability on the transport proper- 
ties, e.g. the magnetoresistance oscillations and match- 
ing anomalies were observed. Further experiments with 
periodic arrays of magnetic dots or holes [18-21] confirm 
that the FM dots create and pin vortices. 

Motivated by current interest and appeal of dealing 
with a new class of physical systems, we consider a peri- 
odic array of FM dots on SC thin film. We assume that 
not only the location of dots is regular, but also their 
magnetic moments are parallel. Such magnetically or- 
dered dot array can be obtained by the field-cooling from 
magnetic Curie temperature to a temperature below the 
superconducting transition temperature T s . It is worth 
to mention that zero-field cooling leads to either random 
dot magnetization for magnetically hard dots, or to an 
antifcrromagnetic ordering if the dots are magnetically 
soft [8,11]. 

In our previous work [16], we demonstrated that a sin- 
gle FM dot with sufficiently strong magnetization and 



large radius generates a vortex or several vortices under 
itself. In a periodic array of dots they must be com- 
pensated by an array of antivortices, since the total flux 
generated by a vortex is zero. Depending on the dot's 
magnetization and size, there might exist several vortex- 
antivortex pairs in the array's elementary cell. For sim- 
plicity, we consider only one pair. Due to the symme- 
try of a square array, the vortices are expected to sit at 
the dot's centers, while the antivortices are located at 
the centers of the elementary cells. However, symmet- 
ric positions may be energetically unfavorable due to the 
inhomogenaous magnetic field distribution over the film. 
Thus, the main problem is to find positions of a vortex 
and an antivortex in an elementary cell. 

In this system, the magnetic field induced by inhomo- 
geneous magnetization of the dot array penetrates into a 
superconductor and generates the SC vortices, whereas 
the magnetic field generated by the supercurrents and SC 
vortices acts on the magnetic system. The method based 
on London-Maxwell equations to study such a system 
was developed in [16]. Under the assumption that the 
sizes of all structures in the problem are larger than the 
coherence length £, the London approximation is valid. 
In section I, we apply the method of the work [16] to pe- 
riodic FM-SC structures. Next, we apply it to a system 
in which a square array of thin circular FM dots is placed 
upon a thin SC film. We conclude with the results and 
discussion. 



I. ENERGY AND MAGNETIC FIELD IN 
PERIODIC HETEROGENEOUS FM-SC SYSTEMS 

A periodic heterogeneous FM-SC system such as 
the regular magnetic dot array upon the SC film 
or a periodic domain structure in a ferromagnet- 
superconductor bilayer can be studied with the method 
described in [16] .The method to study the ferromagnet- 
superconductor bilayer in the continuous limit, where the 
magnetic system's size is larger than effective penetra- 
tion depth, is given in [22]. In this section, we introduce 
a more general method which works well both in contin- 
uous and discreet limits. In doing so, we assume that 
the periodic structures of interest are made of very thin 
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magnetic textures with the magnetization perpendicular 
to the plane and SC films. For simplicity, we consider the 
limit of zero thickness for films with non-zero 2d magne- 
tization and density of superconducting electrons. Their 
energy is calculated over the surface of the SC film. We 
start with the energy of such 2d systems: 



(1) 



where U vv is the vortex energy, U mv is the energy of 
magnetization-vortex interaction, and U m m is the mag- 
netic energy. The terms in (1) are as follows (see [16] for 
details): 
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where eo = </>o/167r 2 A, (j) — hc/2e is the quantum flux, 
A = X 2 L /d s is the effective penetration depth [23], (V<^) 
is the phase gradient induced by the vortex lattice, a and 
b are vector potential and magnetic field at z = 0, re- 
spectively, m is the 2d magnetization. The superscripts 
to and v indicate the contributions of magnetization and 
vortices, respectively. In periodic system all these val- 
ues are periodic functions of coordinates and can be ex- 
panded into Fourier-series. For any periodic 2d function 
f(r), the Fourier expansion is given by 



f(r) = £f G e' 
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where G are the reciprocal vectors of the periodic system 
and A is the unit cell area. Representing all the values 
a, b, m and V</? as Fourier-series (5), plugging them into 
the integral representation of the energies (2,3,4) and us- 
ing equality j e *( G + G )' r d 2 r = A5g.-G', we obtain the 
following expression for the energy per unit cell: 
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The total vector-potential obeys the London-Pearl 
equation 

V 2 A = \aS(z) - 4ttV x (mS(z)) - ^-{V<p)6(z). (9) 
A 2tt\ 

To find the vector-potentials generated by magnetization 
and vortices, we employ the Fourier-expansions of a, 



(V(p) and m. In this article, we are interested only in 
the FM subsystems with magnetization in the z direc- 
tion. The calculations are therefore done along with the 
Fourier coefficient of m z . The Fourier coefficients of the 
phase gradient are given by (see, for example [23]) 
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where Fg = Yli riie tG ' Vi is the structure factor of vor- 
tices, rii and Vi indicate the vorticity and the position 
of the 7-th vortex respectively. In doing so, we get the 
Fourier coefficients of the vector- potentials generated by 
magnetization in the z direction and vortices: 
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Fourier-image of the magnetic field at the interface is 
bo = iG x ac- Substituting the Fourier coefficients of 
the vector-potentials, magnetic fields, the phase gradient 
and the magnetization into equations (6,7,8), we find the 
energy of the periodic structure: 
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(13) 
(14) 
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These equations have rather general character. The only 
assumptions are that magnetic and superconducting sub- 
systems are 2-dimensional and periodic. 



II. THE SQUARE ARRAY OF FM DOTS 

In this section we consider a square array of the FM 
dots on a SC film. We assume that criterion of the vor- 
tex appearance under a single magnetic dot [16] is satis- 
fied. Then the vortices appear under the FM dots, while 
antivortices appear outside them. The existence of an- 
tivortices is obvious since the energy of a system contain- 
ing only vortices would grow as the cube of the system's 
linear size, whereas, in the "neutral" vortex-antivortex 
system the energy grows as the square of its size. We 
consider a simple case in which each dot creates only one 
vortex under them and one antivortex outside (see Fig.l). 
In this section we find the condition for the appearance 
of the vortex-antivortex state and their locations mini- 
mizing the energy. In order to solve these problems, we 
follow the theoretical procedure outlined in the previous 
section. 
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FIG. 1. The top view of the square magnetic dot array is 
as above. The vortices are confined within the dot's region, 
while the antivortices appear outside. 



The geometry of the FM dots is assumed to be the 
same as that in the work [16], namely the magnetization 
of each dot points in positive z-direction, they are circu- 
lar with the radius of each dot equal to R. The lattice 
constant of the square array is denoted as L. The SC 
film and the FM dots are placed at the heights z = 
and z — d > 0, respectively. Later, we transfer to the 
limit d = to simplify the results. For a square lattice 
with the unit cell containing a vortex-antivortex pair, the 
form factor is: 



Fg = cxp (iG ■ ri) — cxp (iG ■ r 2 ). 



(16) 



The set of the reciprocal lattice vectors for the square 
array is given by equation: G = (27m/ L, 2irs/L), where 
n, s are integers; ri = x\x + yiy and r 2 = x 2 x + y 2 y are 
the position vectors of the vortex and the antivortex, re- 
spectively. In order to find the energy of the system, we 
need the Fourier-coefficients of the magnetization. Em- 
ploying the definition of Fourier coefficients for periodic 



systems in (5), the magnetization of a square array of 
FM dots reads 



m zG 
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(17) 



In this result the distance between the dots and SC film 
d is taken to be zero. 

We compute the vortex energy first. Substituting (16) 
into Eq.(13) gives the energy per unit cell: 
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The cosine term in (18) gives the vortex interaction en- 
ergy; whereas, the first term is the vortex self energy. 
In order to extract the logarithmic contribution of the 
vortex self energy, we split the first term in two parts as 
(4tt£o/£ 4 ) EgIVG 2 - V(G 2 (1 + 2AG))]. The first term 
is logarithmically divergent and equals (2eo/^ 2 ) ln(L/£)- 
The second term is convergent and will be left in series 
form. Thus, the vortex energy is found as: 

u m = jr[ ln ^ - ^fv( x / L ) - ^fw(X/L,r 1 /L,r 2 /L)]. 

(19) 

where the functions f v (X/L) f vv (X/L,ri/L,r 2 /L) are de- 
fined by series: 
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, (20) 



where the notation ^ indicates that the term with 
n = s = must be omitted. 

The magnetization-vortex interaction energy can be 
obtained by inserting (17) into equation (14). The re- 
sult is: 
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(22) 
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Now we are in position to calculate the equilibrium 
locations of the vortices and the antivortices. For this 
purpose the total energy of the system should be min- 
imized with respect to ri and r2. For this minimiza- 
tion the magnetization-magnetization interaction can be 
omitted since it does not depend on the vortex coordi- 
nates. Therefore, only the sum 

U — U vv T Umv (^^) 

must be minimized. 



*y 




FIG. 2. The test dot with the vortex and the antivortex. 



For numerical analysis, it is convenient to rewrite the en- 
ergy in terms of dimensionless parameters. To this end, 
we introduce dimensionless variables, namely A = X/L , 
R = R/L, ?i = ri/L and ?2 = v^/L. We also renormal- 
ize the total energy per cell u by dividing it by e v /L 2 . 
We place the center of the test dot at the origin, and the 
positions of the vortices and the antivortices are calcu- 
lated with respect to that point. Let (x\, y\) and (x2, jb) 
be the respective rescaled coordinates of the vortex and 
the antivortex relative to the center of the FM dot. 

With these modifications, we minimize (23) numeri- 
cally with respect to fi and ?2 for different values of 
m(j)o/£ v , R and A. In our numerical calculations, we take 
A/£ = 50. In order for the vortices and the antivortices to 
appear, the effective energy (see Eq.(23)) must become 
negative. Using this condition, we find the phase dia- 
grams for the vortex states, which depend on R/L, R/X 
and m(f)o/e v . The ratio 

e v ~ 9 n s (T)d s In k(T) 1 ' 

is the relative strength of the SC/FM and controlled by 
temperature. In Eq.(24), n m is the density of magnetic 



atoms, d m and d s are the respective thicknesses of mag- 
netic and SC films. S is the value of an elementary 
spin in the magnet and g is the Lande factor and equal 
to 2. n s (T) = n s (T = 0)(1 - T 2 /T 2 ) is the density 
of super-carriers and n(T) = k(T = 0)/y/l -T 2 /T 2 , 
where k(T = 0) = X 2 L (0)/d s £(T = 0). Ln/c(T) varies in 
the range 3 — 6. For typical values of n s (0) — n m — 
10 22 cm~ 3 , d m — 2d s , S = 2 and near the SC transition 
temperature T/T s ~ 0.95, the ratio in (24) can be around 
30. The curves in Fig. 3 indicate when a vortex and an 
antivortex appear for different R/L ratios. Each curve 
separates the regions with and without vortex-antivortex 
pairs. Namely, the region below the curve for particu- 
lar R/L value represents the region without vortex-and 
antivortex pair. Above the curve, there exists vortex- 
antivortex pairs. As seen in the figure, the larger the 
R/L value, the bigger the m4>o/e v ratio. In addition, 
to create vortex-antivortex pair spontaneously, the min- 
imum value of m<j>o/e v increases as R/X ratio decreases 
for each particular ratio of R/L. 




2 - 
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FIG. 3. The phase diagrams of the vortex and the antivor- 
tex state for different R/L. 

Our simulation shows that when the vortices and an- 
tivortices appear on the lines which pass through the 
centers of the dots and make angles (3 and a with the 
horizontal axis respectively. (3 varies between 37r/4 and 
7r, while a varies between 7r/4 and 7r/3, depending on the 
size of the dot, the period of the array and the magne- 
tization. It turns out that a gets closer to 7r/3 as R/L 
increases, whereas fi approaches n. As the dots are placed 
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further away from one another, the vortex moves closer given in the Table I. 
to the antivortex and the dot's center. Some results are 



TABLE I. The position of vortices and the antivortices for different values of the R/X and m0o/e«- The three columns on 
the left are input. 



d / r 
K/L 


it/A 


m(po/s v 


fi 


?2 


P 


a 


0.35 


4 


9 


0.29 


0.41 


180" 


58 


0.30 


4 


9 


0.24 


0.36 


174" 


56 


0.25 


4 


9 


0.20 


0.31 


166 u 


53 


0.20 


4 


9 


0.14 


0.25 


161° 


51 u 


0.15 


4 


9 


0.06 


0.19 


157° 


49 u 


0.35 


2 


9 


0.29 


0.41 


178 u 


58 u 


0.30 


2 


9 


0.24 


0.36 


171° 


56 u 


0.25 


2 


9 


0.19 


0.31 


164" 


53 u 


0.20 


2 


9 


0.12 


0.25 


159° 


52 u 


0.15 


2 


9 


0.04 


0.20 


148" 


48 u 


0.35 


1 


15 


0.28 


0.41 


178 u 


58 u 


0.30 


1 


15 


0.23 


0.36 


171 u 


56 u 


0.25 


1 


15 


0.17 


0.31 


164 u 


53 u 


0.20 


1 


15 


0.10 


0.26 


157° 


50° 


0.15 


1 


15 


0.03 


0.20 


149 u 


48 u 



As seen from our results, the symmetry of the ground 
state is lower than that of the Hamiltonian. The latter 
has point group and it has two sets of rotation axis 
positions: at the centers of the dots and at the centers 
of squares formed by the dots. The residual symmetry of 
the ground state is Z2, with reflection in all specific lines 
of squares. This peculiar symmetry violation stems from 
a highly inhomogeneous magnetic field distribution: the 
screened magnetic field is logarithmically divergent near 
the boundaries of the FM dots. As a result, the vortices 
and the antivortices are dragged into these regions. 

III. CONCLUSIONS 

We considered a periodic array of FM dots on a SC 
film. In the first section, we presented the formalism 
which is modified for the case of periodic structures. In 
this formalism, the problem is formulated as a variational 
principle which allows us to calculate directly the posi- 
tions of vortices. We applied this formalism to the case of 
the square array of the FM dots. We have demonstrated 
that the dots can generate the vortices under themselves 
and antivortices at interstitial locations. Thus, the sym- 
metry of the dot array is spontaneously violated by vor- 
tices and antivortices. For different values of the size of 
the dot and the magnetization, we calculated the posi- 
tions of the vortices and the antivortices and the phase 
transition curves which separate the regions with and and 
without vortex-antivortex pairs. The transition is of the 
first order in our approximation. It turns out that, when 



the ratio R/L increases, then the minimal value of the 
ratio m4>o/s v necessary for the spontaneous creation of 
vortex-antivortex pairs decreases. We assumed that only 
one vortex-antivortex pair per elementary cell of the ar- 
ray is generated. The problem of many pairs has not yet 
been solved. 
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